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0.1 Introduction 

In this note we describe the Chameleon software we developed for the event reconstruction 
of KM3 detector. This software package's developement started as a standalone application 
before the endorcement from the KMSNeT consortium of the SeaTray software framework, 
but it was adapted to it on the course. 

Chapter 1 outlines the techniques we developed for the pattern recognition and the track 
fitting. In Chapter 2, we demonstrate the performance of the Chameleon Reconstruction. 



Chapter 1 

Pattern Recognition and Fitting 



There are two main parts in the reconstruction package: The patern recognition and the 
fitting. The first part is designed for a km^ multi PMT Optical Modules detector (see e.g. 
[2]). It consists of a algorithms for the selection and grouping of hits in track canidates. 

The second is a generic y^ minimizer. This is generic enough to allow for unbiased com- 
parison between different detector designs. 

Before the description of the the pattern recognition and fitting algorithms we give a short 
description of the specific data sets used for this report. 

1.1 Data Sets 



The demonstration of the reconstruction is done on a sample of 6,000,000 MC neutrinos 
produced through the de facto standard tools provided with the seatray, namely nugen for 
the production and g4sim for the simulation. In this data set, no light scattering in water 
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Figure 1.1: The detector geometry used in this study. 



was simulated. ^^K noise was simulated at a rate of 8KHz per PMT, which is more than the 
commonly accepted rate for these tubes of 5.2KHz, since more noise would be a stringent test 
of the filtering out of the noise. For a detailed description of the whole simulation, see [llj. 

The detector used is a small, dense detector of 154 strings (fig. |1.1| ), with string distance 
92 m. Each string has an active length of 570 m., and 20 Multi PMT Optical Modules at 30 
m. distance between each other. Each OM is equipped with 31 3'^ PMTs. 



The produced ^/-spectrum, which corresponds to an £" ^ flux, is presented in fig. L2, while 
the respective quantities for charged current produced muons in fig. |1.3[ 
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Figure 1.2: Generated u spectra. 



A measure of the total number of produced photons is the number of OMs that fired, which 
is shown in fig. 1.5a , for CC events and in |1.5b for NC events. CC events are defined by the 
existence of a primary muon in the particle tree, while NC events by the absence thereof. 

The flat input neutrino spectrum is transformed into the relevant fluxes with the use of the 
"NeutrinoFlux" module (implemeneted originaly for IceCube). Specifically we use functions 
Convent ionalNeutrinoFlux ( "bartol_numu" ) and 
PromptNeutrinoFlux("naumov_rqpm_niimu"). The first simulates the conventional muon flux 
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Figure 1.3: Generated CC fi spectra. 
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Figure 1.4: Angle distributions for /i's.: Black all /i's, Red reconstructed /i's 
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(b) NC events 
Figure 1.5: Numbers of OMs fired. 



(leptons from tt^'s and Kaons), while the second simulates leptons from charm decays. The 
atmospheric flux is the sum of these two functions. 



1.2 Fitting 

An integral part of the pattern recognition are the methods developped for the fitter, so a 
description of the fitting algorithm will necessarily precede the description of the hit selection 
and final reconstruction. 

The fitter is an implementation of a x^ minimization with the help of the C++ Minuit2 
package (embedded within ROOT). The minimization is done according to the geometrical 
model of fig. |1.6[ 




Figure 1.6: The geometrical model used for the minimization. The Cerenkov index used is for 
A = 460nm (43.5°). 

Specifically, if a muon with (pseudo) vertex at R and momentum direction p, emits a 
photon that hit a PMT at x, the total time between the vertex and the hit is given by 



c t^^'P = sin^eV (x - ^)^ - [p • (x - R)]^ + p • (x - R) 
where 9c the Cerenkov angle. The x^ function that is minimized is 



_^exp _ _^data 
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(1.1) 



(1.2) 



where t^^^ is the expected arrival time of the zth photon of eq. (1.1), assuming that the 
measured pulse is the PMT response to Cerenkov light originating from a muon track, tf^^^ 
is the actual measured time, and a^ is the error ascossiated with the zth hit. This error is set 
equal to 2 ns. 

Within the project, two distinct strategies for the estimation of track parameters were 
implemented. They are based on an assesment of the hit residuals wrt the track, as it is 
produced by the fit. 

• Fit with Rejection: This is a recursive fit, which calculates the hit residuals and rejects 
those hits that are Na away from the nominal value of 2 ns. 

• Constant deweight: All hits are deweighted and their a is set equal to a constant param- 
eter. 

The fitter is able to calculate the track parameters in either cartesian coordinates (track 
vertex R(x, y, z) and zenith and azimuth angles) or in spherical coordinates (vertex R(r, 0, (j)) ) 
for a track candidate with at least 6 PMTs hit. 

1.3 Pattern Recognition 

1.3.1 ^^K filter and hit selection 

One of the problems with a sea z^-telescope is ^^K noise. The first part of chameleon needs to 
clean up as much noise as possible, since a x^ reconstruction is very sensitive to noise. 

To use the potential of multi-PMT OMs, the algorithm uses a slightly modified wrt the 
proposed "write up string designs" [2j scheme. 

The first modification is that for each PMT only the first photon is considered to be a hit, 
while the charge is set equal to the total number of photons this particular PMT registered (i.e. 
one has only one hit per PMT). This was deemed necessary until the final time over threshold 
mechanism in the MC is stabilized. The second modification is that not only adjacent PMTs 
are taken into account (see below). 

The ^^K filter is based on a two-layered algorithm, the trigger algorithm and the hit 
selection. 

The trigger algorithm is based roughly on the LI trigger (see [2J), as broadened after a 
MC study of timings for photons. Specifically, for all the hits that are registered on an OM, 
the first thing to do is to select the data set that will be sent to the reconstruction. Since 
signal photons are concentrated around the track that created them while noise is randomly 



produced, a first filter can be based on time coincidences of hits. In fig. (1.7) the differences 
in time of arrival between succesive hits for OMs with 2 or more hits are plotted for noise 
production rate of 8 kHz per PMT. It should be noticed here that although the expected noise 
level for the Multi-PMT design is at 5.2 kHz per PMT, a higher level of noise was used for 
reasons, in order to be on the safe side. If both hits are photons, they are registered in blue. 
It can be seen clearly that when two succesive hits on the same OM are registered within 6 ns 
or less, then they are most probably both real photons. 

The detailed description of these algorithms is as follows: 




(a) Time differences between succesive hits per OM. 

Percentage of noise (orange) and hits (blue) 
1 




(b) Ratios of Cerenkov photons and noise photons wrt the 
total number, for the same sample. 

Figure 1.7: Time differences between pairs of succesive hits per OM. Blue = both hits are photons, 
orange = at least one hit is noise. 
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Figure 1.8: For each OM with more than one photon, choose the 2 photons which are closest in time 
and plot their time difference. One such pair per OM can be plotted. 

• Only use the first photon for each PMT within an event. This is done in order to keep the 
results totally unbiased as to the exact electronics of the system. After the finalization of 
the T/threshold mechanism, the information from the rest of the photons will be used. 

• Start by looping over all the OMs of the event. At this layer consider only OMs with 2 
or more detected hits. 

— After ordering the hits of each OM in time, find which two of them are closest. 
Check their time difference At: If the following condition does not hold 



|At| < 6 ns, 



(1.3) 



the whole OM is discarded, as it is assumed that since the smallest time difference 
is larger than the OM's diameter the photons are too sparse to be real photons, see 



e.g. [6^ p. 60] and fig. (1.7). 

— In case eq. ( |1.3[ ) is satisfied, then perform a local loop over the photons of this 
particular OM, starting from the minimum time pair. Keep only those photons 
that are no further away than 7 ns from the previous photon. 

• After this first filtering, find those OMs that have registered the largest numbers of hits 
(number of PMTs times charge). Only keep the event if the "largest" OM has registered 
hits in at least 3 of its PMTs. These "large" OMs are subsequently grouped by causal 
connection. Each of these clusters of "large" OMs will be used as a basis for track 
reconstruction in the next steps, leading possibly to multiple tracks. 

• Loop over these OM's closest neighbours, i.e. those OMs that are located within 180 
m (^ 3 absorbtion lengths) from the OM with the largest hit. The closest neighbors 
are OMs that did not pass the trigger, so by checking causal connections with the large 
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(a) Time residual for hits with charge equal to 1. 
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(b) Time residual for hits with charge equal to 3. 
Figure 1.9: Time residuals. 
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• 



OMs, more hits can be retrieved. For string distances larger than 180 m this means that 
this search is reduced to same-string OMs. 

Obviously, the number of active (i.e. that have at least one hit on one of their PMTs) 
neighbouring OMs is smaller than the number of neighbouring OMs. The photons on 
the closest neighbours' PMTs are required to obey the condition 

-To< \At\-- <to, (1.4) 

c 

where r is the distance between OMs under consideration, Tq and to are arbitrary con- 
stants and the time difference At is calculated between the current photon's time and 
the time of the first photon of the "large" OM under consideration (For a justification 
see e.g. [7, fig 5.2].) 



In fig. 1.10, the quantity \At\ — - of eq. (1.4) is plotted for a data sample of a small 



detector with string distances at only 90m. As can be seen from the first of these figures, 
for a neighbourhood of i? = 120 m around the largest OM, one can take into account 
all the photons below to, since there the ^^Ks remain a small percentage of the total 
number of photons. Unfortunately, as the radius increases, the local time window, equal 
to the time it takes a muon to travel one such radius ^ R/c, is large enough for more 
"^^K photons to be considered: R/c = 400, 600, 667 ns, for the 3 radii. In the following 
tests To was set to 300 ns for i? = 180 and 200, while it is not important for 120 m where 
the filter is more effective (in the latter case of course there are fewer photons so that 
the next step of reconstruction is harder) . 

The combination of these 2 filters eliminates most of "^^K photons, as can be seen in fig. 



• 



1.11, A major problem whcih should be considered is the fact that as the number of 
OMs under consideration rises, so does the relative number of "^^K hits which survive the 
causality filter. 

What we are left with at this stage is a series of clusters of "large" OMs, together with 
isolated hits from their neighbouring OMs, which are all causally connected. Every such 
group of causally connected OMs is then passed through a fit with rejection method (see 



Q.2) with a N parameter equal to 5 a. The resulting track is used to search again for 
isolated hits. The search is performed on all hits that were not used in the previous 
steps and which are within 2 absorption lengths from the track, on the basis of the time 
residuals of these hits wrt the track. After the search, and if additional hits were found, 
a final fit with rejection at 5 a is performed again. 

1.3.2 Track Parameter Estimation 

One of the main problems with track reconstruction is the fact that (depending on the track 
energy) most of the photons produced by it do not come directly from the muon, but are 
produced through brehmsstralung processes away from it. The result is that the errors ascos- 
siated with the reconstruction are not expected to be gaussian in form, or equivalently, the 
X^ -probability distribution is not expected to be flat. The decision was taken to always pass 
the reconstructed track from a final constant deweight minimization, which aggravates the 
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(b) Photons in a sphere of 180 m. 
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(c) Photons in a sphere of 200 m. 



Figure 1.10: A plot of eq. (1.4). Only photons within a sphere of 120, 180 and 200 m respectively 



are taken into account. Black: survived real photons. Red: Survived "^^K. Blue: Rejected real photons. 
Green: Rejected ^^K. 
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Figure 1.11: The results of the ^^K filtering-out process over all OMs with 3 or more hits plus the 
largest hit of each event and its closest neighbours. String distance 120 m. 
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X^ -probability distribution (by assigning a y^ probability close to 1 or for most of them), 
in order to obtain a more representative pull distribution for the track parameter errors. For 



details, see ^2.1 



1.3.3 Final Track Selection 

As already noted, the above procedure allows for multiple tracks to be reconstructed. In order 
to eliminate ghost tracks, a search is performed amongst the reconstructed tracks for common 
hits. In case such common hits are identified, the track with the smaller overall number of 
hits is deleted if more than 10% of its hits were shared with the track with the larger number 
of hits. 

1.4 The User Point of View 

Thanks to the modular construction of SeaTray, the user only has to decide about the basic 
parameters of the reconstruction and implement them in the driving chameleon module in a 
simple python script. 
The parameters are: 

• The names of input and output files (both *.i3 and *.root files are supported for the 
output) 

• The names of the various containers for particles and photons (input MC particles, 
output reco particles, hits, pulses etc.) 

• The "HitsWeightingMethod" parameter, which specifies the way hits used are going 
to be reweighted during reconstruction. This parameter accepts one of three possible 
values: 

— "constant": This is the default. In this case the error is set by hand in the 
parameter "deweightval" (see below) for all the hits of the track . 

— "residuals": The error is set by the relative time residuals of the hit vs the 
recotrack. 

— "deweight": The deweight and fit method, a recursive method which calculates the 
hit residuals and deweights those hits that are No away from the nominal value of 
2 ns. N is set by the "deweightval" parameter (see below). 

• The "deweightval" parameter, which is used with the "constant" and "deweight" 
fitting methods above. 

The "deweightval" and "residuals" methods were developped for the study of the MC 
data, but were not used in the results shown in this report. 
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1.4.1 Analysis 

The analysis of the output data can be performed with the help of a series of modules that 
export the relevant quantities in chains of root files (see directory AnalysisUtility). A series 
of scripts were written in order to make plots etc (within the scripts directory). 

The implementation of the boost libraries within SeaTray allows for the use of both ROOT 
and python (i.e.. pyROOT) scripts for the same task (see e.g. in the scripts directory files 
timeDiff spyROOT.py and analysis/src/analysis_tree.C). 

In cases where there is need for more computational speed, a relevant C++ class can 
be found in AnalysisUtility, which is also implemented in a way that allows for its us- 
age from within both python and CINT (ROOT's C++ interpreter). For an example see 
Ef f AreaPlotsMaker and its C++ avatar useMakeEf f AreaPlotsClassO, in analysis/src- 
/analysis_tree.C and in pythonic form in angularPlot spyROOT.py. 
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Chapter 2 

Reconstruction Performance 



2.1 Reconstruction Performance 



The efficiency of the reconstruction algorithm is shown in fig. |2.1[ The orange hne represents 
the energy spectrum of ah MC muons that gave hits on at least 6 OMs (Cerenkov photons 
originating from the muon). That is the theoretical minimum of needed OMs for a reconstruc- 
tion. 
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(a) The efficiency of the reconstruction vs the en- 
ergy of CC muons 
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(b) Normalized efficiency of the reconstruction. 
Divided by the number of generated CC muons 
with more than 6 OMs hit in the sample. The 
fluctuation for E < lOOGeV is due to small statis- 
tics. 



Figure 2.1: The efficiency of the reconstruction vs the energy of CC muons. Orange line: Generated 
CC muons with more than 6 OMs hit in the sample. Black line: Generated CC muons with more 
than 10 OMs hit in the sample. Blue line: Contained generated CC muons with more than 10 OMs 
hit. Dashed Black line: Reconstructed muons with more than 6 OMs in the track. Magenta line: 
Reconstructed muons with more than 10 OMs in the track. Red line: Reconstructed muons with 
more than 10 OMs in the track, reconstructed within less than 5° of the MC direction. Green line: 
Reconstructed muons with more than 10 OMs in the track, reconstructed within less than 0.4° of the 
MC direction. 
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This spectrum is compared to the spectrum of reconstructed muons (e.g. the green hne 
represents those tracks which contain at least 10 OMs and their direction is within 0.4° of the 
MC direction). Out of a sample of 5,800,000 neutrinos, 492,147 were CC events with 10 or 
more hit OMs, 253,386 of which were containecQ There were 122,302 reconstructed events 
with 10 or more OMs in the track, 118,930 of which originated from a CC muon. 
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(a) Reconstructed muons with 6 OMs. 



(b) Reconstructed muons with 8 OMs. 
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(c) Reconstructed muons with 10 OMS. 
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(d) Reconstructed muons with 12 OMS. 



Figure 2.2: The efficiency of the reconstruction vs the energy: Reconstructed tracks after cuts, com- 
pared with the respective numbers of "fakes" , i.e. reconstructed muons with a reconstructed direction 
more than 5° off of the MC muon (solid black lines). All plots normalized to the number of Generated 
CC muons with 6 or more OMs hit in the sample. 



Additionaly, fakes are shown in fig. |2.2[ By "fakes" here we designate those reconstructed 
tracks whose direction is 5° or more away from the MC direction. The differences in recon- 



structed and MC 6 and angles are shown in fig. 2.3 



The cut of 10 OMs for the acceptance of a track as well reconstructed follows from plots 
2.2 since it is shown there that a resonable compromise between the reconstruction efficiency 



and fakes is achieved for this cut. This is further corroborated by fig. 2.4, where the difference 



^We define a track to be contained if it passes through the detector, i.e. if the Une of the track crosses at 
least one of the 8 planes that define geometrically the detector. 
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in the reconstructed and the MC neutrinos and muons is shown, before and after the cut. The 
overall median angle for these plots is: 



• Between primary and reconstructed, before cuts, 4.55°. 

• Between primary and reconstructed, after cut, 0.12°. 



• Between MC muon and reconstructed, before cuts, 4.53°, and 

• between MC muon and reconstructed, after cut, 0.11°. 

The differential median angle between MC muon and reconstructed per energy bin is shown 



in fig. 2.5 
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Figure 2.3: Differences in angles between MC muon and reconstructed muon in deg. Black line: all 
reconstructed muons, Red line: muons with at least 10 OMs, Blue line: CC only events with at least 
10 OMs. 
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2.1.1 A note on the asymmetry of fig. |2.3a| and on ghost solutions 



There is an obvious asymmetry in fig. |2.3a between up and down. This is shown separately 



for upcoming (fig. |2.6a ), and downgoing (fig. |2.6b ) tracks. This asymmetry can be possibly 
traced for its most part to hadronic processes, which emit undirectional spherical waves of 
photons. As a demonstration, in fig. |2.7[ the 9 differences are plotted again, but for the red 
line (tracks with at least 10 OMs included), we did not include those tracks that more than 



half of their hits are of hadronic origin. Comparison with fig. |2.3a| shows immediately that 
the number of misreconstructed tracks has dropped by almost one order of magnitude. 

The physics of the asymmetry might possibly be attributed to the detector asymmetry: 
OMs are not spherical symmetric (there is no top PMT). Since then OMs are slightly up-down 
asymmetric (one more PMT in the south hemisphere), the tracks' 9 angles (we remind that 
9 — {) corresponds to a track coming from the nadir) tend to be slightly pulled towards the 
zenith when symmetric light is emmited (as is the case for hadronic processes): the OMs accept 
more light on their south hemispheres from a spherical wave than on their north hemispheres, 
thus interpreting the spherical wave as an upcoming track. Work is in progress to implement 
a hadronic processes module, able to identify and reconstruct correctly these tracks. 



The problem of ghost solutions is of different origin. As it can be seen in fig. 2.8, where 
the same artificial cut on hits of hadronic origin was applied, although the tails are minimized 
after the cut, the kink at around 90° remains. The percentage of the misreconstructed tracks 
around the kink (from 60° to 100°) is 1.04% of the total reconstructed tracks with at least 10 
OMs. 



2.1.2 Quahty of the reconstruction 

The quality of the reconstruction can be assessed by the usual criteria of the pull distribution 
and the x^ probability distribution. The latter is presented in fig. 2^, while the pulls for 9 and 
(j) in fig. |2.10[ These plots were produced with the reconstruction using the Constant deweight 
method (see ^1.2), with a deweight parameter equal to 6 (see below). For completeness, we 
include here the a values for the gaussian fits of these distributions: 



distribution (fig. |2.10a): 



When the distribution without cuts (fig. |2.10a[ left) is fitted, a = 1.549. If only 
the peak of the distribution is fitted (from —2 to 2), then a = 1.097. 



When the (/) distribution with a cut at 10 OMs (fig. 2.10a, right) is fitted, a = 1.162. 
If only the peak of the distribution is fitted (from —2 to 2), then a = 1.024. 



9 distribution (fig. 2.10b): 



When the 9 distribution without cuts (fig. |2.10b[ left) is fitted, a = 1.615. a. If 
only the peak of the distribution is fitted (from —2 to 2), then a = 1.161. 



When the 9 distribution with a cut at 10 OMs (fig. |2.10b| right) is fitted, a = 1.209. 
If only the peak of the distribution is fitted (from —2 to 2), then a = 1.066. 
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(a) Before cuts. 
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(b) After cuts (10 OMs in reco track). 
Figure 2.4: Angle between reconstructed track and muon (left), track and neutrino (right). 
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Figure 2.5: Median between reconstructed and MC muon as a function of logiQ{E). 
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Figure 2.6: 6 difference, downgoing and upcoming. Color code: see fig. 
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Figure 2.7: Demonstration of the up-down asymmetry (see ^2.1.1). Black: MC data. Red: Recon- 
struction, 10 OMs. Up: Declinations. 
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Figure 2.8: Angle between reconstructed track with at least 10 OMs hit and muon (black). The same 
for tracks with less than half of their hits of hadronic origin (red). 
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The naturally ascociated error with each hit (a^ in eq. 1.2) is initially set equal with the time 
uncertainty of 2 ns (more or less equal with the time resolution available). For these plots 
the hits were reweighted with the rather large value of 6 ns, and the justification can be seen 
in plots |2.lT] and 2.12 , which were produced with an error value of 2 ns. The sharp peak at 
of the x^ plot 



2.11 



shows that we underestimate the error. This fact is also reflected to 

This 



the pull distributions (flgs. 2.12) whose a is about 3 or 4, as opposed to flgs. 2.10 



in turn, is a reflection of the highly non-gaussian nature of the errors ascociated with the 
problem, due mainly to the multitude of different interactions involved. For example, the 
number of photons emitted from a track and the form of their distribution in space and time, 
when brehmsstralung and hadronic proceses are involved, is drastically different than the one 
expected from the simple formula used, eq. |1.1| which assumes that all hits are prompt 



photons. 

Therefore, the choice was made to use a larger than expected error for the fltting process, 
since then the pulls' a becomes 1. The apphcation of a reasonable cut (such as 10 OMs), 
diminishes the tails also, which means that in this case the errors on the track parameters as 



calulated by the flt are more representative of the true situation (see flg. 2.10). This, in short, 
is the reason that a large choice for the fltting error is more adequate. 
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Figure 2.9: y^ prob distribution. Black: All tracks. Red: reconstructed with 10 or more OMs. 



The neutrino effective area vs energy is plotted in flg. 2.13, where we compare the effective 
area for all reconstructed tracks (blue line), to the effective area when only those tracks with 
at least 10 OMs are taken into account (black line). As a measure of the real efficiency of the 
reconstruction, we also plot (the red line) the effective area without the fakes (i.e. those tracks 
that are reconstructed more that 5° away from the true direction). 
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Figure 2.10: Pulls for (/), 6>. Left: all tracks, Right: tracks with 10 OMs 



25 



1 chi2prob | 


chiprob 




Entries 134786 
Mean 0.1461 
RMS 0.2816 
Underflow 
Overflow 


10= 

10^ 
10= 


- , ^-u. < 1 < < < 1 < < < 1 < < < 1 < 


1 











0.2 


0.4 


0.6 


0.8 


1 







Figure 2.11: x^ prob distribution for a constant deweight fitting with error 
Black: All tracks. Red: reconstructed with 10 or more OMs. 



2 (see text, §. 2.1.2) 
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Figure 2.12: Pulls for O^cf) for a constant deweight fitting with error = 2 (see text, §. 2.1.2). Left: all 
tracks, Right: tracks with 10 OMs 
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Figure 2.13: u effective areas. Blue: all reconstructed events. Black: only events with at least 10 
OMs. Red: events with at least 10 OMs, reconstructed within 5° from the MC direction. 
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Figure 2.14: Same as 2.13, but only for upgoing muons {z > 90° 
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Chapter 3 

Depth Studies 



The measure of the quahty of a detector is its "sensitivity" [13], and here we expore its relation 
to depth. The procedure for the calculation of the sensitivity is based on the Feldman-Cousins 
technique [12j. 

The estimation of the atmospheric v backgrounds was based on the MC files, described 
already in §1.1[ The atmospheric [x backgrounds on the other hand, pose a difficult problem, 
due mainly to their large numbers - and the resulting need for resources. Within the present 
work, a rough MC estimation was done for them, while the more quantitative full Monte Carlo 
study will follow soon. 

After a short description of the method we use for calculating the fluxes in the flrst section 



of this chapter, the results are presented in §3.2 

3.1 Fluxes 

The flrst step in calculating the point source sensitivity, is the evaluation of atmospheric and 
astrophysical v fluxes. The outline of the calculation is as follows: 

• Extract the events' zenith and azimuth. Calculate declinations and right ascensions (see 



^3.2 for a discussion of some relevant technical details and flg. 3.4 for a sky map). 

• Extract the weights, as calculated from the simulation (for an explanation see |10|). 

• Calculate the signal rate R\ as 

n = ^t^T, (3.1) 

where N the total number of events, O^ the OneWeight parameter, T the live time in 
seconds and 

- '■''">" (3.2) 





^i - 




E^ 


is the signal flux. 








• Calculate the relevant background rate 


as 






R\ 
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(3.3) 




Figure 3.1: Expected Atm. neutrino rate (as per eq. 3.3). The nature of the peaks at ±1 can be, at 
least partiahy, attributed to ghosts, see §2.1.1[ Black: All reconstructed events. Blue: Reconstructed 
with at least 6 OMs. Red: Reconstructed with at least 10 OMs. 



where now $J is a bartol flux [9J (see ^1.1 and fig. 3.1). The factor of 2 is due to the 
counting of both us and i>s. 

All the above factors are related to the irreducible background of atmospheric neutrinos, 
which is practically unchanged with depth. The only factor that depends on depth is the flux 
of atmospheric muons. Due to limitations of time, the full Monte Carlo study of atmospheric 
muon fluxes will be performed at a later stage. 

From the user point of view, the task is divided in two parts. The first part is implemented 
in the Ef f AreaPlots module, which outputs the data needed for the calculation to a chain of 
ROOT files. The second part of the calculation is finished by a series of scripts manipulating the 
data (see esp. Ef f AreaPlotsMaker and the functions therein, and the relevant implementation 
scripts, e.g. analysis_plots.C and eff AreapyROOT.py) 



3.2 Sensitivity 

We calculate the point source sensitivity by using the straightforward fixed bin method, where 
we use an angular bin of the order of the angular resolution (For a basic definition of terms see 
[8j). The main idea is that signal events will cluster around the direction of the astrophysical 
neutrino source under consideration (within of course the angular resolution of the detector), 
therefore producing an excess of events over the uniform background. 

The background is made of two parts: the atmospheric neutrino part and the atmospheric 
muon part. The former is already included from the MC source files that were reconstructed. 
For the latter, a full reconstruction is not easy to be performed because of time and CPU 
constraints (work is underway in this direction). Until the inclusion of the full Monte Carlo 
data for the atmospheric muon fluxes, a general argument can be used: On one hand, the 
atmospheric muon intensity, as a function of sea depth, is a well documented quantity (see e.g. 
P,J, (fig. 1-3) and discussion thereafter, or [3J, fig. 15). On the other, the number of expected 
atmospheric muons is already given from the existing simulations. From a preliminary study 
of atmospheric muons at 3500 m, using the program mupage, the estimated muons that reach 
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the detector would be of the order of 8 x 10^ events per year. This number is corroborated by 
(fig. 1-3) of [Ij, for a detector of approximate surface Ikm? . 

The detector produces fakes (here the word "fakes" signifies downgoing events that were 
misreconstructed as upcoming) at a known percentage of input events. In fig. 3.3a the ratio 



of misreconstructed events over the total produced downgoing MC neutrino events is shown. 
This ratio can be used for an evaluation of the number of misreconstructed atmospheric muons, 
since their energy spectrum, as a function of depth, is known. The result of this exercise is 



shown in fig. 3.3b 




Figure 3.2: Energy spectrum of atm /i's. Red: Depth = 5000m. Black: Depth = 4500 m. Blue: 
Depth=3500 m. 
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(a) Ratio of misreconstructed "fakes" over to- (b) Expected atm muon "fakes" over one year 
tal number of MC events. period. 

Figure 3.3: Black: All reconstructed events. Blue: Reconstructed with at least 6 OMs. Red: Recon- 
structed with at least 10 OMs. 

The total number of estimated fakes, at 10 OMs cut, is 94754, for a live time of one year. 
For a uniform distribution of the fakes over the whole lower part of the sky and an angular 
bin of radius 0.4°, this is reduced to 2.3 fakes per bin per year. It should be noted that, 
due to low statistics, t he nu mber of fakes is overestimated slightly (observe the fluctuations 
around IQ^GeV in fig. 3.3b), and averaged over the whole sky, something that is expected to 



be amended in the near future, after a complete reconstruction of the full data set. 
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As a comparison to the above number, the mean number of expected atmospheric v back- 
ground events per angular bin of 0.4° per year is 4.4. 

These numbers are in good general agreement with the relevant MC simulations published 



so far: see e.g. fig. |3.6[ where the atmospheric muon intensity at 3500 m is of the same order 
of magnitude as the neutrino background. As a qualitative rule then, one can extrapolate 
the information about the atmospheric background from its two sources, calculating the total 
background (Atm /i's + Atm z/'s), as follows: If at 3500 m, the //-background is expected to 
be on the average half the v background, then at 2500 m, their ratio is expected to be roughly 
2, and conversely at 4500 it will be 0.11 and at 5200 0.04. 



The plot for the sensitivity for a cone of 0.4° around the MC direction is shown in fig. 3.5 



A calculation of the sensitivity of our detector for various //-background scenarios and for 2 



different cuts on the number of OMs is presented in Table [3T| for a declination of 53°, near 
the minimum. The left column is the ratio of the atm-/x background over the v background, 
i.e. when this quantity is 1, there are equal numbers of v and [i background events. 



Table 3.1: Sensitivities (GeV" 



X 10 ^^), with an estimated atm /i background. 
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Figure 3.4: A sky map of reconstructed muons in equitorial coordinates declination 5 and right 
ascension r". Upper: locally upcoming muons {z > 90°). Middle: downgoing muons {0° < z < 90°). 
Lower: Total. Not a very useful plot, but it does display nice colors. 
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Sensitivity for MC v *s with z«n = 90 to 180 10 r«co OMs cut 




Figure 3.5: One year sensitivity at the Nestor site, for a cut at 10 OMs. Green line: no atm 
background. Red Line: An estimate for atm mu background equal to atm v background. Blue line: 



An atm mu background equal to 10 times the atm u background. For a discussion see section 3.2 
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